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A variety of neutrino flavour conversion phenomena occur in core-collapse supernova, due to the 
large neutrino density close to the neutrinosphere, and the importance of the neutrino-neutrino 
interaction. Three different regimes have been identified so far, usually called the synchronization, 
the bipolar oscillations and the spectral split. Using the formalism of polarization vectors, within 
two-flavours, we focus on the spectral split phenomenon and we show for the first time that the 
physical mechanism underlying the neutrino spectral split is a magnetic resonance phenomenon. In 
particular, we show that the precession frequencies fulfill the magnetic resonance conditions. Our 
numerical calculations show that the neutrino energies and the location at which the resonance takes 
place in the supernova coincide well with the neutrino energies at which a spectral swap occurs. The 
corresponding adiabaticity parameters present spikes at the resonance location. 

PACS numbers: 14.60.Pq„97.60.Bw 



I. INTRODUCTION 

Intense theoretical activity has unravelled the complex- 
ity of neutrino flavour conversion in media. The neu- 
trino propagation in a star like our sun is at present 
well understood as a resonant flavour conversion phe- 
nomenon due to the coupling of neutrinos with ordi- 
nary matter, the Mikheev-Smirnov-Wolfenstein (MSW) 
effect [TJ [2]. This confirmation has in particular been 
brought by the SNO [3J, KamLAND g] and recently by 
the Borexino [5] experiments. Neutrino oscillations in 
massive stars turn out to be a non-linear phenomenon 
because of the importance of the neutrino-neutrino inter- 
action, as first pointed out in [51 [7] and also later on [5]- 
110] , In particular, collective flavour conversion phenom- 
ena have been shown to emerge |llj , the understanding of 
which has already required a large theoretical effort (see 
e.g. 0ESHHII3II3II5H13) and [Mir27] for a review). 
Besides, in such explosive environments the presence of 
shock waves [23j [28H33] and of turbulence [34H37] also 
induces new flavour conversion effects. The inclusion of 
all these features has significantly increased the complex- 
ity in the computation of neutrino propagation in such 
media. Fully understanding the flavour conversion phe- 
nomena and their interplay with key unknown neutrino 
properties, such as the neutrino mass hierarchy, the third 
unknown mixing angle and leptonic CP violation [HO [H] 
is one of the major goals in this domain. While impor- 
tant progress has been achieved, several aspects still need 
further clarification. 

Currently, three collective flavour conversion regimes 
induced by the interaction of neutrinos with themselves 
are identified, usually called the synchronization regime, 
the region of bipolar oscillations and the spectral split 
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phenomenon. The first occurs very close to the neu- 
trinosphere : any flavour conversion is frozen due to 
the neutrino-neutrino interaction |llj . Then the vac- 
uum term triggers an instability [HI Q21 EH [25] (depend- 
ing on the neutrino luminosities and the hierarchy) that 
produces flavour oscillations. It appears to be associ- 
ated with a rapid growth of the derivative of the mat- 
ter phase introduced by the neutrino-neutrino interac- 
tion [35]. The onset of such an instability is different 
in calculations based on the single-angle approximation 
versus a multi-angle treatment [24) . Bipolar oscillations 
have been shown to correspond to electron anti-neutrino 
and neutrino pairs converting to pairs of neutrinos and 
anti-neutrinos of other flavours |13j . Finally, as first in- 
dentified in [11] the neutrinos undergo a complete flavour 
conversion depending on their energy. This implies the 
swap of the electron neutrino spectra with the ones of 
muon and tau neutrinos, below or above a critical value 
called the split energy. The spectral split has been fur- 
ther studied in [TTJ W\ (Figure [XJ> and Refs. [HI [22] 
have shown that such a phenomenon is even more com- 
plex since multiple spectral splits show up, depending 
on the neutrino flux ratios at the neutrinosphere and on 
the core-collapse supernova explosion phase. While sev- 
eral aspects have now been identified, a comprehensive 
understanding of the spectral split phenomenon is still 
missing. 

Extensively used in the literature, the formalism of 
neutrino polarization vectors constitutes a good tool to 
gain a better understanding of neutrino flavour conver- 
sion phenomena. A long time ago it was pointed out 
that neutrino oscillation in vacuum could be pictured 
as the precession of a neutrino flavour polarization vec- 
tor around an effective magnetic field, depending on the 
neutrino masses and mixings [39]. Later on [40] ST] 
the case of the resonant (adiabatic and non-adiabatic) 
flavour conversion in matter has been discussed, in con- 
nection with the MSW effect and the solar neutrino 
deficit problem. Recently such a formalism has been 
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used in the presence of the neutrino-neutrino interaction 
Hamiltonian in supernovae first in and subsequently 

in nana ca EDi 03]. 

Using such a formalism in Ref. |42j the synchronization 
regime has been investigated in the context of the Early 
Universe. The bipolar oscillation regime has been put in 
relation with a flavour pendulum in [111 115] and with a 
gyroscopic pendulum in [T3]. In [13113] the authors make 
the hypothesis that at the end of the synchronization and 
of the bipolar regimes, the neutrino evolution follows a 
collective precession mode until neutrino densities are low 
and at the late stage of this precession solution the step- 
wise swapping occurs. In [20j . an explicit adiabatic solu- 
tion of the two flavour neutrino evolution equations with 
the neutrino-neutrino interaction is built in the comov- 
ing frame. The existence of such a frame had already 
been pointed out in the early work [TTJ . It is shown that 
analytical solutions exist and present a behaviour like 
the one of the neutrino spectral split. Approximate con- 
served quantities, in particular the total neutrino lepton 
number for the two-flavour case, have been identified and 
have been shown to explain the split energy observed in 
the neutrino fluxes [2D] . This idea has been further inves- 
tigated in [17] . Besides the synchronization and the pure 
precession modes already mentioned, a self-induced para- 
metric resonance mode is discussed in |18j . Many of these 
studies have brought very valuable insights into the col- 
lective effects engendered by the neutrino self-interaction 
and have given a vision of what actually is seen in the 
full numerical calculations. Note that, in order to allow 
for these useful analytical solutions, approximations are 
often made : (i) neglecting the matter term; (ii) using a 
constant coefficient for the neutrino self-interaction term 
or with a simplified time dependence; or (iii) employ- 
ing a single neutrino or anti-neutrino energy or neutrino 
box spectra. An application of the polarization vector 
formalism to the three-flavour case, requiring the SU (3) 
basis, is performed in [12] and it has been shown that 
two conserved quantities also appear in this case [43l [44] . 
Finally, adiabaticity of the neutrino evolution has been 
particularly discussed in [501 HH1 HI] ■ 

In this work we approach the problem of the neutrino 
spectral split using the matter basis and the polariza- 
tion vector formalism within two neutrino flavours. We 
give the explicit analytical expressions for the polariza- 
tion vectors and the effective magnetic field, that in our 
case includes both the neutrino-neutrino interaction, the 
vacuum contribution and also the matter term. In the 
matter basis we follow the full neutrino evolution in the 
region where the neutrino-neutrino interaction effects are 
dominating, and focus, in particular, on the spectral split 
phenomenon. We show, for the first time, that a mag- 
netic resonance mechanism is the underlying process gov- 
erning the flipping of the polarization vectors and there- 
fore the spectral split. We point out that the flipping - 
giving rise to a spectral swapping - occurs at the neutrino 
energies for which the precession frequencies meet the 
magnetic resonance condition. We show the behaviour 
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FIG. 1. Neutrino fluxes as a function of neutrino energy at 200 
km near the neutrinosphere, in the case of inverted hierarchy 
and within two neutrino flavours. The spectral split due to 
the neutrino-neutrino interaction is apparent as a swap of the 
electron neutrino with the v x neutrino spectra. The results 
correspond to neutrinos (left figure), where the split energy 
is at 7.6 MeV, and anti-neutrinos (right figure), with a split 
energy at about 2.4 MeV. 



of the associated generalised adiabaticity parameters. 

The paper is structured as follows. Section II re- 
calls the polarization vector formalism for two-neutrino 
flavours in connection with the neutrino Hamiltonian de- 
scribing their propagation in a core-collapse supernova. 
Section III gives the equations to follow the evolution of 
the effective magnetic field and of the polarization vec- 
tor. Section IV establishes the link between the spectral 
split phenomenon and a magnetic resonance. It presents 
the numerical results for the neutrino probabilities and 
adiabaticity parameters in the region where the neutrino- 
neutrino interaction effects are important, and shows 
that the resonance condition is indeed satisfied when the 
spectral split occurs. Section V is a conclusion. In the 
Appendix, expressions to calculate the second derivative 
of the neutrino-neutrino interaction Hamiltonian in the 
multi-angle treatment, are provided. 



II. THE POLARIZATION VECTOR 
FORMALISM AND NEUTRINO EVOLUTION IN 
A SUPERNOVA IN THE MATTER BASIS 

Let us start by reminding the well known fact that the 
quantum evolution of a system described by an Hermitian 
2x2 matrix 



H = 



H21 H22 



(1) 



with H 2 \ — H^ 2 i can be associated with an effective mag- 
netic field 




2TZ(H 12 ) 
-2X{H 12 ) 
(i?H — H22) 



(2) 



with 1Z and I indicating the real and imaginary contribu- 
tions of the off-diagonal matrix element. The evolution 
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of the system can be seen as the precession of an effec- 
tive spin around B. Indeed, if we define the polarization 
vector P 



(3) 



with a the Pauli matrices, the evolution equatiorj^] can 
be written as 



rfP 

dx 



= B x P 



(4) 



with precession frequency u = |B|. 

Let us now consider the specific case of the neutrino 
evolution in a core-collapse supernova. The correspond- 
ing Hamiltonian in the flavour basis is : 



(5) 



where the non-diagonal vacuum contribution is UKU^ 
with K = diag(k 1 , k 2 ), U is the Maki-Nakagawa-Sakata- 
Pontecorvo (MNSP) unitary matrix, relating the flavour 
to the mass basis [46]. The other two contributions cor- 
respond to the matter ternj^] Hm ^ — VZGfPc, coming 
from charged-current interactions of neutrinos with the 
electrons in the medium, and the neutrino-neutrino in- 
teraction term Hli [14]. We will follow the neutrino 
evolution in the matter basia^J as done in Ref.[5S]. In 
such a basis the Hamiltonian is : 



K = U^H U ">U = diag(k 1 ,k 2 ) 
with the two-neutrino flavour Hamiltonian 



(6) 



(7) 



and U the MNSP matrix depending on the matter angles 
and phases, (k±,k2) the matter eigenvalues. The evolu- 
tion equation reads 



i— = Hip 
dx 



(8) 



with ip = (vi,v 2 ) the matter eigenstates and the matter 
Hamiltonian involving the derivatives of the U matrix 



dx 



More explicitly one gets 



H 



h + Qi 



(Sk 12 +8Q 12 ) 



_. (Sk 12 +SQ 12 ) p* 
' 2rr 1 



(9) 



(10) 



1 Note that in the manuscript we will use x to denote time. 

2 Note that the neutral current contribution is taken out since it 
only contributes with a term proportional to the unity matrix 
and is therefore irrelevant for neutrino oscillations. 

3 From now on, the symbol O will always indicate that the quantity 
(in this case O) is calculated in the matter basis. 



with Ski2 — k\ — k 2 the diagonal contributions coming 
from the U derivative 



with i = 1,2, and we define 




(11) 



— &i 2 I U ^ 



(12) 



The off-diagonal terms are given by the generalized adi- 
abaticity parameters [25) 



2 7re"* 12 



(Sk 12 + SQ 12 
2ire l5112 




(13) 



Sk 12 (6k 12 + 5Q 12 ) V dx / 12 
We take the most general expression for U : 



U 



1 

e l ? 



cos 9 siaO \ ( e~ lai 
-sine? cost? Mo e 



(14) 



where j3 and cti,a 2 are the Dirac and Majorana phases 
in matter respectively. The neutrino evolution equation 
becomes explicitly [25] : 



dx V v 2 



+ psm 

sin20 



sin20 



+ it 



(3cos 




where we use the dot to indicate d/dx. The Majorana 
phases (ai,a 2 ) are set to zero without loss of generality 
since they do not influence the neutrino flavour conver- 
sion. 

We now apply the polarization vector formalism to 
Eq.(10). Following Eqs.(I][2), one gets the effective mag- 
netic field corresponding to the matter Hamiltonian 



B 



-^.x(r) 



k 



12 



-^sin26> 

26 I (K)i 

Sk 12 - /3cos20 , 



The magnetic field x- and y-components only depend 
upon the Hamiltonian off-diagonal element, while the z- 
componcnt is the difference of the diagonal ones. From 
Eq.([3]) one can construct a neutrino polarization vector 
in the matter basis 



— 21(^1^2 ) 



(17) 



"2 



whose z-component tells us the v\ or v 2 content of the 
neutrino wavefunction. The x- and y-components are as 
usual non-zero in presence of mixings. 
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Now solving Eq.(15l is equivalent to determining the 



evolution of the polarization vector Eq.( 171 under the ac- 
tion of the effective magnetic field Eq.(16), i.e. Eq.Q. 
In this work, we retain all contributions to the Hamilto- 
nian including the matter term. In order to follow the P 
evolution we need to calculate the matter angle derivative 



(18) 



H ee - H^j Sku - (H ee - 5k 12 

4|ff eM |5fci2 

and the phase derivative : 



1 



l(H efl )TZ(H efl ) - K(H efl )l(H^) (19) 



Eqs.( 18|19 l involve the derivatives of the diagonal H e 



H^^ and off-diagonal H e „ Hamiltonian matrix elements 
in the flavour basis Eq.([7]) (see [S3] for details). Note 
that, if matter only is included, only the derivative of 
the matter angle is relevant, while the presence of com- 
plex off-diagonal contributions in the Hamiltonian in the 
flavour basis introduces also the j3 phase and its deriva- 
tive. 



III. THE EVOLUTION OF B AND P 

In order to establish a connection between the spec- 
tral split phenomenon and the magnetic resonance phe- 
nomenon let us now define relevant quantities. First of 
all the total precession frequency of the polarization vec- 
tor Eq.(17) around the effective magnetic field Eq.(16) 



Up 



IBI 



(20) 



while the angle defining the evolution of B with respect 
to the (XOY) plane is defined as 



tan 77 



5k 



12 



B 



2iJi 



(21) 



In order to follow the evolution of the system it is inter- 
esting to define the rate of change of the B direction with 
respect to such a plane, namely 



u)„ = rj = 2/„ 



Sk 12 \H 12 \ - 2Sk 12 



K{H 12 )K(H 12 ) 



-1(H 12 )1(H 12 ] 



(22) 



with the normalization factor /„ = [4\Hi 



One sees that to follow Lb v both the derivative of 5k\ 2 
and of r are needed. For the former, one gets 



Sk 12 = B z = Sk 12 - ft cos 29 + 2(59 sin 29 



(23) 



with 



Sk 12 = 



(u*hWu) - (frirWifj (24) 



while one gets for the second derivative of $ 



(25) 



with = \H elx 



n{H eil )Tl{H ell ) +l{H eil )l{H eil ) 



The other important quantity to 



determine for Eq. ( 22 1 is 

5ki 2 8k 12 r = — ( 5k 12 8k 12 + 5ki 2 8k 12 ) T 



7r(u^H if) U+ U^H( f) U,fPtJ ^(26) 



which depends upon the second derivative of Eq.([5j). 
Such a derivative involves two contributions, since the 
vacuum term is constant 



(27) 



The first term depends on the explicit matter density 
profile used, that we take as a power-law. Its contribu- 
tion to Eq.(27) is straightforward. As far as the neutrino- 



neutrino interaction term is concerned, the corresponding 
second derivative depends on the use of a single-angle ap- 
proximation versus a multi-angle treatment. Here we give 
its expression in the former case since it is the approxi- 
mation we employ in the following. The relations valid 
in the multi-angle case are given in the Appendix. The 
neutrino-neutrino interaction Hamiltonian in the single- 
angle approximation^] is given by 



H, 



(/) _ 



F(x)G(p) 



(28) 



where p indicates the density matrix for two neutrino 
flavours and the geometrical factor is 



F(x) 

If = V2G F /(2nRl) and 



[!-<?(*)] 



g{x) = \i- 



(29) 



(30) 



and R v the radius of the neutrinosphere. The expression 
for the non-linear contribution is explicitly 



G(p) = 



a=e,fi 



(31) 

where L Va is the neutrino flux at the neutrinosphere for 
a neutrino "born" as an a flavor. The second derivative 
of H vu includes contributions from both the derivative of 
the geometrical factor and of the density matrices, i.e. 

H vv = F(x)G{p) + 2F{x)G{p) + F(x)G{p) (32) 



4 Note that in our calculations the neutrinos are emitted with the 
same angle taken to be 0° with respect to the neutrinosphere as 
done e.g. in Ref. |14| . 
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To calculate £j ri we first need the obvious second deriva- 
tive of the geometrical term 



F(x) 



\g{x)- 2 + (1 - g(x))g(x)- 
+ 3x 2 R- 2 (l-g(x))g(xy 



(33) 



The second contribution to Eq. ( 32 1 can be calculated as 
follows 

(34) 



where 



F(x) = -f F x 6 Ri 



3r>2 



(35) 



Finally the third one can be evaluated using 



IV. THE CONNECTION BETWEEN THE 
MAGNETIC RESONANCE AND THE 
SPECTRAL SPLIT PHENOMENON 

In its simplest realization, the magnetic resonance phe- 
nomenon consists in a spin-flip occurring in presence of 
a constant and a time varying magnetic field. Let us 
consider a spin precessing around a constant magnetic 
field B = i?o z i which can be taken along the z-axis, 
with a precession frequency given by ujq. The spin is 
also subject to a second magnetic field Bi, located in 
the (xOy) plane, which sinusoidally varies in time with w 
frequency. Such a magnetic field acts as a perturbation 
that can however significantly impact the dynamics of the 
system depending on the relation between Aw = w — ujq 
and Wi, the precession frequency around Bi. If Aw ^ 
the time varying B x has little effect on the spin evolu- 
tion. However, if the resonance condition Aw ~ is met 
and Aw <C wi, the presence of Bi can fully flip the spin. 
The amplitude for the resonant flip conversion follows a 
Breit-Wigner distribution w 2 /(w 2 + Aw 2 ) which is maxi- 
mal when Aw ~ 0, with the resonance width being given 
by u}\. On the other hand, if Aw is a multiple of wi the 
system is off resonance while if the frequency difference 
is a fraction of it the system is close to resonance and the 
spin inversion is partial [45] . 

Let us now define three relevant frequencies that are 
useful to establish our connection between the neutrino 
spectral split phenomenon and the magnetic resonance 
phenomenon (Figure [2]). The effective magnetic field as- 
sociated with the Hamiltonian Eq. ^ in the flavour basi^] 




FIG. 2. Analogy between the spectral split phenomenon and 
the magnetic resonance phenomenon, based on the polariza- 
tion vector formalism. The figure depicts the effective mag- 
netic field B, its z-component Bz and the component lying 
in the (XOY) plane Bj_ with the corresponding frequencies 
(see text). 



has components in the (xOy) plane proportional to the 
off-diagonal H ell matrix element which depends on the 
vacuum and neutrino-neutrino contributions; while the 
third component is generated by the matter term (that 
gives a large contribution) as well. 

What we would like to argue now is, that the inversion 
of the neutrino polarization vectors is due to the fact that 
the precession frequencies meet the magnetic resonance 
criteria. To this aim, we switch to the matter basis to fol- 
low the system evolution. Using the polarization vector 
formalism (Section II) we identify the Bz = £>o compo- 
nent Eq.(16l for the B magnetic field. We also identify 



5 Note that we use small x,y,z (capital X,Y,Z) letters when in 
the flavour (matter) basis. 



the Bi time varying magnetic field with B^. We now first 
show that different flavour conversion regimes are well 
accounted for, if one follows the average of the magnetic 
field instead of the field itself. Let us have a closer look 
at the time dependence of the magnetic field direction. 
Figure [3] plots the fast variations of the r\ angle Eq.(21) 
showing that B undergoes fast oscillations with respect 
to the (XOY) plane. The numerical results presented in 
this work have been obtained by using the vacuum oscil- 
lation parameters \Sm 2 \ = 2.4 x 10~ 3 eV 2 , 8 — 9° and the 
matter density profile pb = 1-5 x 10 8 (10/x), with x be- 
ing here the distance in the supernova, in units of g.cm 
and of km (for x). The neutrinosphere R v is taken at 10 
km and the corresponding neutrino fluxes are assumed to 
be of Fermi-Dirac type with average energies (E Ve ) = 12 
MeV, (E„ K ) = 15 MeV and (E Ux ) = 18 MeV. Equipar- 
tition of energy is assumed, with a total luminosity of 
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4 x 10 51 erg.s" 1 , While we present one set of our numeri- 
cal results, we stress that we have performed calculations 
for an ensemble of initial conditions, varying the neutrino 
average energies and luminosities at the neutrinosphere. 
The quality of the results we obtain is the same for all 
initial conditions considered, so that we show here one 
case only. 

Our calculations of the ratio of the ui^ Eq. ( [22] ) and lip 
Eq.(20l frequencies, i.e. Q v /Qp shows that such ratio is 
always much larger than 1 : the magnetic field direction 
changes rapidly compared to the precession of the neu- 
trino polarization vector around it. Therefore let us see 
how well the system is described if we replace B by (B) . 
To do such a comparison, in order to obtain our (B), 
we perform a numerical average of the magnetic field B 
over 10 km. If we define the a angle describing the Bj_ 
evolution in the (XOY) plane, 



tan a 



(By) K av (Y) 



(B X ) Iav(T) 



(37) 



where the subscript av means that we are taking the 
average of the corresponding quantity lZ a v(Y) = (TZ(Y)) 
I(r)) and the associated angular velocity, 



andl a „(r) 



1 



u)„ — a = 



r I 2 

1 av I 



(Y)1l av (Y) - K av (Y)X av (Y) 



The uj a frequency involves the derivative of the Y factor 
Eq.( |26[ ). In our numerical calculations we obtain that 
the TZ av (Y), lZ av (Y) and therefore the Q„ frequency are 
very close to zero. The second quantity of interest is the 
precession frequency around Bz, namely u) = \(Bz)\ = 

\5ki2 — /3cos20|. Another relevant precession frequency is 
the one around the Bi, namely 



|B3°| 



B 



x 



((/3sin2<?) 2 ) + (40 2 



(39) 

Since the contribution coming from (A9 2 ) comes out to 
be small, the quantity o>i is essentially determined by the 
phase derivative. Finally one can also define the same 
quantities as in Eqs.(20) and (22) but this time for the 



average magnetic field, namely 



Wp 



(B) 



(40) 



The rate of change of the angle defining the direction of 
(B) with respect to the (XOY) plane is now given by 



tan 77 



(B; 



(5k 



12/ 



Bf 



(41) 



with Bf 



((B x ), (By))- We show in Figure|3|an exam- 
ple of the evolution o f th e i] angle defining the evolution 
of (B) given by Eq.(41| for a 5 MeV neutrino (black 




FIG. 3. Time evolution of the r\ angle Eq.(| 21[) de fining the di- 
rection of the effective magnetic field B EqTIl6l with respect 
to the (XOY) plane (grey curve). The effective magnetic field 
is associated with the full matter Hamiltonian Eq.(10l in- 



cluding the vacuum contribution, the matter term and the 
neutrino-neutrino interaction (in the single-angle approxima- 
tion). The fast oscillations are driven by the derivative of the 
P phase. The black curve shows the angle with respect to the 
(XOY) plane giving the evolution of (B) instead of B Eq. |4l| 
The results are obtained for a 5 MeV neutrino and inverted 
hierarchy. 



We have therefore numerically verified, that the evo- 
lution of our system is well described if one uses the 
(B) = ((Bx), (By), (Bz)) instead of the magnetic field 
itself. To this aim we have performed two calculations for 
the neutrino evolution and determined the corresponding 
oscillation probabilities. In the first calculation we have 
solved the precession equation Eq.Q with B. In the sec- 
ond we have determined (B) making a numerical average 
over 10 km (as mentioned above) and then solved Eq.Q 
using the average of the effective magnetic field instead 
of the field itself. Figure [4] presents the comparison of 
the survival P(vi — > Vi) probabilities with i = 1, 2 for the 
matter eigenstates, obtained when the neutrino evolution 
is determined by either B or (B). The system evolves 
near the neutrinosphere where the neutrino-neutrino in- 
teraction is important, with the evolution determined by 
the full matter Hamiltonian Eq.(10). Results are shown 



curve), instead of B (grey curve). 



both for neutrinos and for anti-neutrinos of two different 
energies. One can see that the synchronization, bipolar 
and spectral split regimes can be easily identified in both 
cases. Note that the different behaviours of the prob- 
abilities shown in Figure [4] engender the characteristic 
spectral swap of the electron with the other v x flavour 
fluxes (Figure []}. One can see that, when averaging the 

fast variations due to (3, the main features of the syn- 
chronized, the bipolar and the spectral split behaviours 
remain unchanged. Note that the two descriptions agree 
qualitatively and quantitatively very well until 120 km 
for the neutrino, and 80 km for the anti-neutrino cases. 
Afterwards the evolution of the systems following the av- 
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erage of B and B itself start presenting some differences. 
However, looking carefully at the survival probabilities 
both in the matter and in the flavour basis (Figure[5]), one 
can see that the region where the two descriptions com- 
pare well is sufficient for our purposes since discrepancies 
appear at the very end of the spectral split behaviour in 
all cases. This convince us that the results we will present 
in the following are not an artefact of following the aver- 
age of the effective magnetic field. Therefore, from now 
on, we will replace B by (B). All the results of the next 
section are obtained following such a prescription. 



Neutrinos 



Anti-neutrinos 




FIG. 4. Evolution of the third component of the polarization 
vector Eq.( 17 1 as a function of distance within a core-collapse 



supernova, determined by using the magnetic field B Eq. (161 
(black curves) or its average (B) (gray curves). Results are 
shown for a 5 and a 10 MeV neutrino (left-figure) and for a 
5 and a 2 MeV antineutrino (right figure). The numerical 



calculations include the full matter Hamiltonian Eq.( 10 1 with 
vacuum mixings, neutrino coupling to matter and to neutri- 
nos (in the single-angle approximation). The neutrino energy 
dependent behaviour characteristic of the spectral split phe- 
nomenon (giving rise to a complete flavour conversion and a 
spectral swap) can be seen. 



V. THE OCCURRENCE OF THE SPECTRAL 
SPLIT AND THE MAGNETIC RESONANCE 
CONDITION 

The main goal here is to show that the magnetic reso- 
nance conditions : 



Auj 



UJq = 

< 1 



CO 



= 
< 1 



(42) 



are verified for the neutrino or anti-neutrino energies and 
at the same location in the supernova for which the spec- 
tral split occurs. We emphasize that verifying that the 
spectral split is a magnetic resonance phenomenon can 
be done both in the frame identified by the flavour basis 
and in the one given by the comoving frame. In the flavor 
basis the effective magnetic field has indeed a small time 
varying component Bj^ in the (xOy) plane; the magnetic 
resonance conditions are those given on the left of (42 1. 



basis and the average of B, (Bj_) appears as static and 
the magnetic resonance conditions are those on the right 
side of ( 42 ) . For convenience we show the fulfillcmcnt of 



the magnetic resonance conditions in the comoving frame 
identified by the average of the magnetic field in the mat- 
ter basis, i.e. (B). 

Let us follow the different phases of the neutrino evolu- 
tion with the polarization vector formalism in the matter 
basis. At early times, near the neutrinosphere, the mat- 
ter term is so large that the polarization vector evolution 
is essentially determined (Bz) ~ (Sku)- At this stage 
it is essentially precessing around this component and 
nothing happens from the point of view of flavour con- 
version. As shown in [25], at the end of the synchroniza- 
tion regime, the onset of bipolar oscillations is triggered 

by a rapid growth of the /3 phase. This implies that the 
Bf components - related to the off-diagonal elements 
of the matter Hamiltonian - start being important, as 
does the phase contribution to (Bz)- hi fact the phase 
growth introduces an oscillating degeneracy between the 
off-diagonal Hamiltonian matrix element and the differ- 
ence of the diagonal matrix elements. As a consequence 
the ratio of such elements determining the generalized 
adiabaticity parameter, become comparable, as discussed 
in [25] . In the present formalism such parameters are 



given by cot7y (see Eq.(41 )) : 



IB1 



\(2H; 



12/ 



(Bz) 



\(Sk 21 +SQ 21 ) 



(43) 



During the bipolar oscillations the polarization vector 
evolution is determined essentially by the two compo- 
nents (Bj, By). However when taking the average, since 

the average is around zero (note that the derivative 
of angle turns out to be much smaller than the phase 
derivative), the effective magnetic field associated with 
the matter Hamiltonian stays in a XZ plane. This is 
the co-rotating frame already discussed in the literature 
QTJ E3] • After some distance the polarization vector 
enters in a new phase of its evolution : the magnetic res- 
onance regime. In this phase, depending on the neutrino 
energy, if the resonance condition Au = uj — ujq w is 
satisfied then a full inversion of the neutrino polarization 
vector occurs. In the comoving frame in the matter ba- 
sis such a condition readf0 Gj ~ and near resonance 
ui/ui\ <C 1. 

Figure [5] presents the evolution, in the flavour basis, of 
survival probabilities for neutrinos and anti-neutrinos in 
inverted hierarchy, in correspondence with the |T| evo- 
lution. The latter quantity, multiplied by the factor 
{8k\2 + 8Q 12) / 5k\2, is also shown to emphasize the role of 
the (SQ12, coming from the U'dU/dx contribution to the 
diagonal elements of the matter Hamiltonian. Results are 



On the other hand, in the frame identified by the matter 



6 Note that in the corotating frame u> = U3 — tJrj' 
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FIG. 5. Evolution of survival probabilities of the flavour 
eigenstates (upper), of the corresponding generalized adia- 
baticity parameters V Eq.(43l (lower) and of the same quan- 
tity multiplied by the factor (<5fci2 + SQi2)/5ki2 (middle fig- 
ures). The results are obtained by solving Eq.( |I5| l for the full 
Hamiltonian comprising the vacuum contribution, the mat- 
ter term and the neutrino-neutrino interaction (in the single- 
angle approximation). The case is inverted hierarchy. In each 
figure the two curves correspond to different neutrino ener- 
gies, namely a 5 MeV (black) and a 10 MeV (gray) neutrino 
or anti-neutrino. 
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FIG. 6. Three dimensional contour plot of the electron neu- 
trino survival probability (in the flavour basis) as a function 
the neutrino energy and distance within the supernova. The 
colors encode the electron neutrino survival probability. The 
yellow-green region corresponds to a survival probability of 
0.5 so that all neutrino energies below the split energy of 7.6 
MeV do not undergo any spectral swap, while all energies 
above the split energy undergo a spectral swap. The black 
curve superposed shows when the magnetic resonance condi- 
tion in the matter basis is met : uj = Eq. ( 42 1 (see Figure 
|8». 



presented for two different neutrino and anti-neutrino en- 
ergies chosen to show the different behaviour of the gen- 
eralized adiabaticity parameters for neutrinos that do not 
(do) undergo the spectral split phenomenon and a swap- 
ping of the fluxes depending on their energy. One can 
see from the figure that the corresponding |T| has a very 
different behaviour in presence of the swaps (i.e. the anti- 
neutrino and the 10 MeV neutrino cases) with respect to 
the absence (the 5 MeV neutrino case) . In particular, for 
the former cases at the location of the resonance spikes 
appear in the adiabaticity parameter evolution, at 80 km 
for the 10 MeV neutrino, and around 60 km for the anti- 
neutrinos. Note that the second spike at 110 (130) km 
for a 5 (10) MeV anti-neutrino is not associated with any 
significant change because the neutrino-neutrino interac- 
tion term has become negligible and the second condition 
Q/Q\ <C 1 is not satisfied (see Figs. [8] and [9]) . 

Figures [6] (neutrino case) and [7] (anti-neutrino case) 
present three dimensional contour plots of the neutrino 
survival probabilities in the flavour basis for different neu- 
trino energies, as a function of distance within the super- 
nova. Figures [8] and [9] show the regions where the reso- 
nance uj (black lines) and near resonance Q/Qi -C 1 
(band within white curves) conditions are met in the mat- 
ter basis. The coloured area outside the white shaded re- 
gion is a region off-resonance oj/uji S> 1. By comparing 
Figure [6] with [8] and Figure [7] with [9] one sees that : i) 
the same neutrino and anti-neutrino energies for which 
the magnetic resonance condition is fulfilled also undergo 



a spectral swap; ii) the inversion of the neutrino flavour 
polarization vector occurs at the same distance in the 
supernova where the resonance occurs. One can also 
see that the resonance region in the anti-neutrino case 
is much smaller than in the neutrino case in agreement 
with the faster decrease of the anti-neutrino survival (or 
appearance) probabilities in the flavour and matter ba- 
sis, compared to the neutrino case (see Figure [5]). Finally 
the upturn of the black and white curves in Figures [8] and 
[9] shows a region where the condition uj m is fulfilled 
while Q/uji ~ 1. No inversion of the neutrino flavour 
polarization vectors is observed in such a case, in corre- 
spondence with the second spike observable in the adia- 
baticity parameters T (Figure |5j. When the precession 
frequencies are way out of resonance and the role of the 
neutrino-neutrino interaction clearly negligible, the sys- 
tem evolution is determined again by the component, 
given by Sk±2- 



VI. CONCLUSIONS 

We have investigated collective flavour conversion ef- 
fects, within two-flavours, induced by the neutrino- 
neutrino interaction in a core-collapse supernova. Us- 
ing the neutrino polarization vector formalism associated 
with the matter Hamiltonian and the neutrino matter ba- 
sis, we have focussed upon the spectral split phenomenon 
and established a connection with a magnetic resonance 
phenomenon. In such a formalism, the third compo- 
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FIG. 7. Three dimensional contour plot of the electron anti- 
neutrino survival probability (in the flavour basis) as a func- 
tion the anti-neutrino energy and distance within the super- 
nova. The colors encode the electron anti-neutrino survival 
probability. The yellow-green region corresponds to a survival 
probability of 0.5 so that all neutrino energies below the split 
energy of 2.3 MeV do not undergo any spectral swap, while 
all energies above the split energy undergo a spectral swap. 
The black curve superposed shows when the magnetic reso- 
nance condition in the matter basis is met : Q — Eq. ( 42 1 
(see Figure |9j. 



nent of the effective magnetic field is given by the dif- 
ference of the diagonal matrix elements of the matter 
Hamiltonian comprising the matter eigenvalues and the 
phase derivative induced by the presence of complex off- 
diagonal terms in the neutrino flavour Hamiltonian. The 
X- and Y-components of the effective magnetic field de- 
pend upon the derivative of the matter phase and angle 
respectively. By defining precession frequencies around 
the average of the Z-component and of these perpendic- 
ular components, our numerical calculations of the po- 
larization vector evolution have allowed the identifica- 
tion of a new phase in the evolution : the magnetic res- 
onant regime. This is encountered when the magnetic 
resonant criteria concerning the precession frequencies is 
met. Our results show this very interesting feature : the 
neutrino energies for which the resonant criteria are ful- 
filled are the same energies undergoing the spectral split 
phenomenon, occurring at the same location in the su- 
pernova. When the system is in such a regime the cor- 
responding generalised adiabaticity parameters present 
spikes, while they present a smooth behaviour for the 
neutrino energies that do not meet the resonant and do 
not undergo any spectral swap. Once the connection be- 
tween the spectral split and a magnetic resonance phe- 
nomenon established, it is clear that the fulfillment of the 
magnetic resonance conditions can be done either in the 
flavour basis, or in the comoving frame (identified here 
by the average of the magnetic field in the matter basis). 
While we have used the latter for convenience, picking 
one or the other is, obviously, just a question of choice. 




FIG. 8. Connection of the spectral split phenomenon with 
the magnetic resonance phenomenon for neutrinos: The col- 
ored blue region shows the neutrino energies and the super- 
nova location for which the resonance conditions Hi « and 
Cj/Cji <C 1 (on resonance) Eq.(|42[) are both fulfilled. The black 
curve presents the location of the magnetic resonance corre- 
sponding to the fulfillment of the ui ~ condition. The white 
line separates the area where we are on resonance from the 
one off-resonance given by ui/uii 3> 1. Note that the quantity 
uj\ gives the magnetic resonance width. The figure shows that 
the resonance criteria follow a line, one around 75 km where 
the flipping of the neutrino flavour polarization vector occurs, 
that then stays flat and then turns upward at around 120 km. 
Note that in this second region the polarization vectors are 
not flipped since the criteria uj/uii <C 1 is not met anymore 
(the neutrino-neutrino interaction has become negligible). 



In this work the numerical results have been obtained 
assuming equipartition of neutrino energies at the neu- 
trinosphcre and inverted hierarchy since for this case the 
neutrino-neutrino interaction gives rise to single splits. 
While the results presented in this work correspond to 
a set of initial conditions in terms of neutrino average 
energies and luminosities at the neutrinosphere, results 
of the same quality are obtained for a large set of initial 
conditions. We think the mechanism we have been iden- 
tifying to be general and valid also in the case of multiple 
splits and of three neutrino flavours. Such investigations 
will be the object of future work. 



VII. APPENDIX 

Here we give expressions for the second derivative of 
the neutrino-neutrino interaction Hamiltonian valid in 
the case of a multi-angle treatment. In this case the 
non-linear Hamiltonian reads 



rrm.il. _ V%Gf ST^ 

vv ~ 2nRl ^ 



dq' d cos 8' (1 — cos 9 cos 9') 



[ Pu J q ',8')L u J q ')-plJq',0')L,J q ')} 



(44) 
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FIG. 9. Connection of the spectral split phenomenon with 
the magnetic resonance phenomenon for neutrinos: The col- 
ored blue region shows the neutrino energies and the super- 
nova location for which the resonance conditions Co w and 
uj/uji <C 1 (on resonance) Eq.( |42[ | are both fulfilled. The black 
curve presents the location of the magnetic resonance corre- 
sponding to the fulfillment of the ui ss condition. The white 
line separates the area where we are on resonance from the 
one off-resonance given by w/wi 2> 1. Note that the quan- 
tity uj\ gives the magnetic resonance width. The figure shows 
that the resonance criteria follow a line, one around 65 km 
where the flipping of the neutrino flavour polarization vector 
occurs, that then stays flat and then turns upward at around 
100 km. Note that in this second region the polarization vec- 
tors are not flipped since the neutrino-neutrino interaction 
has become negligible. 



with 



COS Qrn.aT — 



^l-{R v /xf 



(45) 



Using the Leibniz Integral Rule : 



OZ Ja{z) Ja(z) ° z 



+f(b(z),z 
-f(a(z),z) 



db{z) 

dz 
da(z) 

dz 



(46) 



one can show that the derivative for Eq. ( 44 ) is [25 



Tjm.a. 



V2G F 
2-kRI 



V/ dq' (R(e',q')-S(x)T(q')) 
„ Jo 



with 

R{9',q') 



S(x) 



[ dcos0' (-«(1-cos0cos0') 

J cos 8 max 

[[H,p V! Xq',e>)\ L^(q') 

+ [H,p^q>,e>)]*LzJ.q')]) (48a) 



x^Jl-{R w /x) 2 
T(q') = {1 - cos9 cos9 max )[p Va (q' ,6 max )L Ua (q') 



(48b) 



-PIM' \°max)L Pa (q')} 



( f) 

The second derivative of Hi, v is given by: 
V2G 



(48c) 



-S(x)Ttf)-S{x)ftf)) (49) 



Its evaluation requires the first derivatives of the R, T 
and S functions. For R we get : 



R(e', q ' 



dcose 1 {R 1 (e',q') + R 2 (d',q'))+R 3 (6',q' 

(50) 



where 

Rl{6',q') = -ifg, {^H^, Pv& {q , ,e , )\ L^{q') 

+ [HV\p 9 Jq',6')]* L Q Jq')) (51a) 

R2(0',q') = -%fv (-1 [H {f \ [H«\p u& (q',e')]] L v Jrf) 

+z [R^\[H^\ P ^q\d')^ L 9s {q')) 

(51b) 

R0,q')=iS{x)fe m ^ ([H( f \p v M\e')\ L„Jq') 

+ [RV\p p Jq',6')]*L Q& {q')) (51c) 



with 



fe> = (1 - cos (9 cos 6>') . 

For the other two one obtains 

3 R^ 



S(x) = Ri 



with 



and 



(47) 



x 4 g(x) x 6 g(x) 3 
g(x) = ^\-{R u /xf 

t{q>) = T x {q>) + T 2 {q>) 



(52) 



(53) 



(54) 



(55) 
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with 



-pi & (q',O max )L D ^(q') (56a) 
r 2 (g') = -*/9 m „. [[# (/ W(g',0 ma:c )] L„ ft (gO 

W>1* W<?')1 (56b) 
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